Weighted Scale-free Networks in Euclidean Space Using Local Selection Rule 
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A spatial scale-free network is introduced and studied whose motivation has been originated in 
the growing Internet as well as the Airport networks. We argue that in these real-world networks a 
new node necessarily selects one of its neighbouring local nodes for connection and is not controlled 
by the preferential attachment as in the Barabasi- Albert (BA) model. This observation has been 
mimicked in our model where the nodes pop-up at randomly located positions in the Euclidean space 
and are connected to one end of the nearest link. In spite of this crucial difference it is observed 
that the leading behaviour of our network is like the BA model. Defining link weight as an algebraic 
power of its Euclidean length, the weight distribution and the non-linear dependence of the nodal 
strength on the degree are analytically calculated. It is claimed that a power law decay of the link 
weights with time ensures such a non-linear behavior. Switching off the Euclidean space from the 
same model yields a much simpler definition of the Barabasi- Albert model where numerical effort 
grows linearly with N. 
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Studying the structure of communication networks is 
important in their own right because it helps in under- 
standing the network, its traffic flow distribution which 
in turn helps making the communication process more ef- 
ficient. Over last several years many real- world networks 
are being studied with much interests. The examples 
vary from social networks, electronic networks and bio- 
logical networks Q, 0, Q • An important subclass of these 
networks are spatial networks, i.e., those embedded in 
the Euclidean space. Two most important examples of 
these networks are the electronic communication network 
like the Internet 0, 0, , which is a transport network 
of electronic data packets as well as the public transport 
system of Airport networks 0, • 

The common property of these two networks is their 
highly inhomogcncous structures. This inhomogeneity is 
reflected in their nodal degree distribution (degree A: of a 
node is the number of links attached to it) which is ob- 
served to follow the power law distribution: P(k) ~ k~ lk . 
It was Barabasi and Albert (BA) who first recognized 
that indeed there are many other real- world examples in 
social and biological networks having similar structures 
0. Since these networks lack a characteristic value for 
the nodal degree, they are called as the Scale-free Net- 
works (SFN) pUjfl. 

BA observed two key characteristics for these net- 
works: (i) they are growing networks (ii) existence of an 
inherent 'rich get richer' mechanism which ensures that 
large degree nodes grow at higher rates. In BA model a 
network grows by addition of new nodes which get con- 
nected to the growing network using a linear attachment 
probability. In this paper we question the general neces- 
sity of this 'linear attachment' rule. We argue that at 
least for spatial networks like the Internet and the Air- 
port Network existence of such a rule for the expanding 
network seems to be highly implausible. 



If you buy a new computer and would like to connect it 
to the Internet what do you do? If it is a home computer 
you use a modem and a dial-up telephone line to connect 
it to the nearest router of the telephone company. If it 
is your office computer, it gets a connection to the office 
router, which is eventually connected to the nearest node 
of the Internet service provider (ISP) your organization 
had opted for. If your office is a part of your university 
or organization, several routers are used and in the au- 
tonomous systems level they are connected to the nearest 
ISP again. Therefore the new nodes of the Internet in any 
level are always connected to the existing local nodes of 
the network. A person in Chile would hardly give any 
extra importance to large ISP hubs in Tokyo, Stuttgart 
or Chicago rather than small providers in his locality. 
Considering the whole world-wide Internet network it is 
apparent that the new nodes pop up randomly in space 
and time without any spatial correlations. 

Similar arguments can be put forward for the airport 
network as well. A new airport in some remote city in 
some country is quite likely to be connected to the neigh- 
boring airport first by direct non-stop flights and very 
unlikely to have cross-country inter-continental flights to 
begin with. Recent studies on IATA airport network data 
also revealed that weights Wij of the links in this net- 
work either defined by the number of passenger seats or 
by the Euclidean distance between successive stops have 
non-trivial variations |lfj| . The nodal strength defined as 
Si = 'EjWij varies as (s(k)) oc k@ with (3 > 1. A number 
of models have been proposed to put forward explana- 
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tions for such non-linearity 

In this paper we study a growing network on Euclidean 
space where new nodes are added one by one and are con- 
nected to the neighboring nodes of the growing network. 
We show that even such a network has scale-free degree 
distribution. In addition, these networks have non-trivial 
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FIG. 1: The picture of a network of N = 2 nodes, each 
node is a randomly positioned point on the unit square and 
is connected randomly to one end of the nearest link. 

dependence of average strength on the nodal degrees. 

A spatial network is grown on a two-dimensional 
space whose nodes are the points at randomly se- 
lected positions within an unit square on the x — 
y plane (with periodic boundary conditions) |l4| . 
Let {(xi,yi), (22,2/2), {xn,Vn)} represent the co- 
ordinates of the of N randomly distributed points within 
this space where each co-ordinate is an i.i.d. number 
drawn within {0, 1}. The growth of the network starts 
with a pair of nodes 1 and 2 connected by a single link. 
Nodes labeled 3 to N are then introduced one by one and 
are connected to the growing network. 

Let at some intermediate stage the network have t links 
and t + 1 nodes. For connecting the (t + 2)-th node to 
the network, the nearest link center is selected. One of 
the two end nodes of the nearest link is then chosen with 
equal probability and is connected to the new node to 
create the (t + l)-th link (Fig. 1). This is executed by 
keeping the local informations into memory. The unit 
square area is divided into a lattice of size y/N x y/N. As 
the network grows the serial numbers of the links whose 
centers are positioned within a lattice cell are stored at 
the associated lattice site. To find out the nearest link 
center one starts from the cell of the (t + 2)-th node and 
searches the lattice cells shell by shell till the nearest link 
center is found out. When t ~ N only the nearest shell 
needs to be searched. 

After the network has grown to N nodes, the degree 
distribution P(k, N) is calculated. From a direct mea- 
surement of the slope of the log P(k, N) vs. logfc plot 
the exponent 7^ is estimated to be 3.00 ± 0.05, which is 
close to BA value. Moreover, a scaling of P(k, N) is also 
studied (Fig. 2(a)): 

P(k,N) cx N-^ik/N^) (1) 
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FIG. 2: (a) Finite-size scaling analysis of the degree distribu- 
tion for three system sizes, N = 2 14 (solid line), 2 16 (dashed 
line), 2 18 (dotted line) and for three values of the number of 
outgoing links m = 1, 2 and 3 (from left to right), (b) Same 
plot as in (a) but scaled with m -1 ' 67 . 

where rj = 3/2 and ( = 1/2 are used to obtain the best 
data collapse giving 7fc = rj/C, = 3. The network so gener- 
ated has a tree structure. However networks having more 
general structures with multiple loops can be generated 
with m links coming out from every incoming node and 
are attached in a similar way to the first m neighbor- 
ing links in increasing order. In Fig. 2(b) we show that 
degree distributions for m = 2 and 3 can be scaled as 
P(k, N)m~ ls7 cx N- r 'g{k/N^). 

The link lengths are also measured and the ir p roba- 
bility distribution is calculated as studied in 14] . Let 
T>(£)d£ denote the probability that a randomly selected 
link has a length between £ and £ + d£. For a given Pois- 
son distribution of N points let us first calculate the first 
neighbor distance distribution. Consider a point P at an 
arbitrary position. The probability that its first neigh- 
bor is situated at a distance within r and r + dr (which 
can be done in TV — 1 different ways) and all other N — 2 
points are at distances greater than r + dr is: 

P(r)dr = {N -I) 2nrdr (1 - Trr 2 )^ 2 (2) 

In the limit of N — > 00 it can be approximated that 
N — 1 m N — 2 Ki N and since the average area per 
point decreases as 1/N, nr 2 is very small compared to 1. 
Therefore (1 — Trr 2 ) N ~ 2 is approximated as exp(— irNr 2 ). 
Therefore in the limit of N — > 00 the the probability 
density distribution is: P(r) ~ 2nNr exp(—irNr 2 ) or in 
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FIG. 3: Finite-size scaling analysis of the link length distri- 
bution D(£,N) for three system sizes, N = 2 14 (solid line), 
2 16 (dashed line), 2 18 (dotted line). Fig. (a) shows the plots 
on a linear scale where as in (b) the double logarithmic scale 
has been used for the same data but binned on exponential 
scale. 

the scaling form: 

P{r)/y/N(x [rVN] exp(-7r[r\/iV] 2 ) (3) 

where the scaling length 1 / y/~N is the linear extent of the 
average area per node. 

In our case the number of nodes N in the system is not 
fixed rather grows with time. Therefore the link lengths 
£ also decrease as time progresses. After some time the 
total collection of links has a mixture of many different 
lengths. Since intially there were only few nodes the early 
links are very large and may be as big as the box size, 
where as the last few links are very small and have lengths 
of the order of £ ~ 1/yN. Effect of the mixing has 
been two fold as observed numerically. The distribution 
of small link lengths (later stage) up to £ m » 2.h£ Q is 
different from Eqn. 3 but still follows a scaling form: 

T>(£, N) oc </Nr(£VN). (4) 

The scaling function fits very well to the generalized 
Gamma distribution: G(x) = Ax a exp(— bx c ) where, 
A w 12.6, a w 1.4, b « 3.11 and c « 0.83 (Fig. 3(a)). 
On the other hand the long early links contribute to a 
power law tail as T>(£) ~ for £ m < £ < y/2/2 with 

lw « 3.0 (Fig. 3(b)). 

Like weighted airport network [lOj the strength s, of a 
node i is measured as the sum of the Euclidean lengths 



FIG. 4: (Color on-line) On an unit square a six link network 
is drawn and centers of the links are marked. The whole 
space is then partitioned into six different Voronoi cells around 
these centers. Cells are marked by dropping 2 16 randomly 
selected points in this space and coloring all points in a cell 
by a particular color which are nearest to the corresponding 
link center. 

of all links meeting at i: Sj = Sj^y, l%j being the length 
of the link between the nodes % and j and a is a continu- 
ously tunable parameter. This parameter generalizes the 
model to take care of situations where the link weight 
may even vary non-linearly with the Euclidean distance. 
E.g., the route that an aircraft takes for flying between 
two airports is quite often greater than the length of the 
geodesic path between them. For example the ASIANA 
airlines flies from Seoul to New Delhi through the air 
space of Bangladesh. Also in Mobile ad hoc network 
(MANET) the power Pj required to maintain the range 
Ri of transmission of each mobile element varies nonlin- 
ear ly with the range: Pj cx Rf, with 1 < a < 6. 

We first observe that given a connected network of 
t links, the whole space is partitioned into t non- 
overlapping Voronoi cells, each cell surrounds the center 
of a link |l5j (Fig. 4). The center of the link is at a 
minimum distance from all points within this cell. The 
probability that a randomly selected point is within a 
particular cell is equal to the area of the cell. Since dif- 
ferent Voronoi cells have different areas the probability 
of selecting a link center is in general non-uniform. For 
a two-dimensional Poisson Voronoi tessellation the cell 
sizes follow a Gamma distribution whose width scales 
as 1/t 16] . Therefore though for finite t the cell sizes 
are non-uniform, for t — > oo the cell size distribution is 
similar to a Delta function implying that all cells have 
uniform size. In this limit the probability that a partic- 
ular node of degree k will be linked is k/2 times the cell 
size - which gives rise to the linear preferential attach- 
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FIG. 5: Verification of the formuias for (a) 7.5(a) = 1 + 4/(2 + 
a) and (b) (3(a) = 1 + a/2 (solid line) with the numerical 
results (circles) for TV = 2 14 . 

merit as in BA model. Therefore for a very large network 
(N — > 00) in our model the node selection probability is 
different from that in BA model at early times (t small) 
where as it asymptotically converges to the BA prefer- 
ential rule in the limit of t — ► 00. We conclude that the 
leading behaviour of our spatial network model is like the 
BA model. We have already seen that even for finite N 
the degree exponent 7^ « 3.0 as in the BA model. 

Let us assume that the area of every Voronoi cell in the 
network with t links is uniform and is equal to X/t. This 
implies that the i-th node with degree ki (t) has probabil- 
ity ki(t)/(2t) to be linked with the new (t + 2)-th node. 
The factor 2 comes from the fact that one of the two 
end nodes of every link is selected with equal probabil- 
ity. Therefore dki(t)/dt = h(t)/(2t) and kilt) = (t/i) 1/2 , 
exactly similar to the BA model resulting P(k) cx fc~ 3 . 

Again, since the area of each cell is 1/t the typical 
length £ij of the t + 1-th link to be connected to the node 
i is proportional to (fcj(i) /t) 1 ^ 2 . Therefore the rate of 
increase of the strength of the i-th node is: 

ds 4 (t)/dt = (dh(t)/dt).(l%) oc (h{t)/t) 1+a / 2 (5) 

which is proportional to (ti)~ 1 / 2 ~ a / i . On integration 
over t from t = i to t: 

»«(') - "«« « - ra/2 )- < 6 > 

The value of Sj(i) is estimated by the average strength of 
the i-th node when it was introduced and connected to 
an arbitrary previous node j,j = 1, i — 1 

*(») cx r 1 / 2 -/ 4 (7) 
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FIG. 6: Variation of the exponent 8(a) with a for N = 2 1 . 

When t is large in Eqn. (4) the term i~ Q / 2 is ignored 

ai(t)«r 1 / a - 0, / 4 [t 1 / 2 -°/ 4 + c] (8) 

Writing i in terms of Sj (t) 

i oc [tV2-«/4 + c ] 4 /( 2+a) h(t)]- 4 /( 2+Q ) (9) 

Since the serial number i and the strength s are both 
nodal quantities they are dependent and their probability 
distributions obey P(s,t)ds = —P(i,t)di, the negative 
sign is because as i increases s decreases. Using P(i, t) = 
1 jt we get the time dependent probability distribution of 
the nodal strengths as: 



(X [^/2-«/4 + c ]4/(2+a) r_^ s -l-4/(2+a)] 

2 + a 

oc f-irks- 1 -^ (10) 

Therefore 7 s (a) = 1 + 4/(2 + a) and using the general 
relation j s = 7*//? + 1 — [8j and using 7^ = 3 we get 

/9(a) = 1 + a/2. (11) 

The exponent of the distribution of Wij — £fj varies as 
r ) w (a) = 1 + 2/ a. A number of checks have been done to 
verify these results. For a large network with N = 2 14 
the strength distribution and its average are calculated 
for 21 different values of a equally spaced between and 
2. Results are found to be quite consistent with the above 
formulas of -f s (a) and (3(a) (Fig. 5). 

Therefore the message is that a nontrivial value of 
(3 > 1 indeed can be obtained in a weighted network 
when the link weights decrease systematically as time 
elapses. In the above model the link weight decreases 
as (it)~ a / 4 . This is consistent with the airport network 
data where early airports in big cities like London, Paris 
etc. have very high strengths. Because of the fact that 
these early airports still survive over more than a century 
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implies that they are connected with strong links of high 
passenger traffics as well as long distance links. There- 
fore if the early links have maximal weights and if they 
decrease inversely with a power of time, that can result 
/3> 1. 

A simple way to check this idea is to study the gen- 
eralized BA model itself, starting from a single link. 
The weights of the links are assigned by hand: A link 
which becomes a part of the network at time t carries a 
weight t~ a . Then if the degree of the largest hub {k max ) 
grows as f (/i = 1/2 for the BA model) then it can 
be shown that: for a < /j, (3(a) — 1 and for a > /x, 
(3(a) = a, I ii. The weight distribution exponent varies 
with a as 7„(a) = 1 + 1/a. We observe that these re- 
sults are very much consistent with the results in [l2| . 
In this model each link comes in with a constant weight 
and at the same time m! randomly selected links gain 
additional weights w at each time step. As a result old 
links increase their weights as time increases. On simu- 
lating the Bianconi model network of size N = 2 14 we 
calculated that the average weight of a link introduced 
at time i decreases as i~ a which is exactly their result 
that the weight of a link increases with time as (t/i) a . 
In comparison to the conclusion of th at a global re- 
organization of weights (as in and yields (3 > 1, 
we claim that ultimately the global reorganization or any 
other dynamics has to ensure a power law decay of the 
link weights like t~ a to achieve (3 > 1. Similar conclusion 
has been drawn in [l2T | . 

To study the correlation between the nodal degrees 
and the link weights Barrat et. al. calculated how 
the product of the degrees kikj of the two end points of 
a link varies with the link weight my 

(w i:i ) cx [kikjfW (12) 

A direct calculation of (wij) averaged over different con- 
figurations generated for different values of a for our spa- 
tial network shows that variation is almost linear (Fig. 
6). 

We further observe that our model can be generalized 
when the new node selects one of the t links randomly 
with a probability where S is a continuously 

tunable parameter. Therefore in the limit of <5 — > — oo we 
retrieve the above model. On the other hand for S = 
the underlying Euclidean space is switched off and the 
networks are generated in a much simpler way: To add 
the (t + l)-th link, a link j is selected out of t links with 
uniform probability l/t. The new (t + 2)-th node is then 
randomly connected to one end of j with probability 1/2. 
Therefore a node of degree k has the probability k/(2t) 
to be linked with the new node which clearly satisfies 
the requirement of linear attachment probability of the 
BA model. Therefore we claim that the algorithm: the 
new node selects one of the links with uniform probability 
and gets connected to its one end with probability 1/2 is 



exactly the BA algorithm. Computationally there is a 
lot of advantage, it takes a CPU linearly proportional 
to N as already observed in 0, IrH . Iisj . The network 
so generated clearly has a tree structure. However, loops 
can also be generated easily by connecting each new node 
with to distinct links, each link is attached with the same 
rule. The resulting network is exactly the BA network 
with m outgoing links from each node. 

Finally it is observed that the case of the non-linear 
attachment rule in BA model when the probability of 
attachment varies as fc e cannot be obtained by general- 
izing this model. Let the two end nodes of the randomly 
selected link have degrees k\ and k 2 then depending on 
their degrees one of them is selected with the following 
probabilities: 

p(k 1 ) = k{ /{k{ + and p(k 2 ) = k\j(k\ + kl) (13) 

with p(k\) + p(k 2 ) = 1. The BA model corresponds to 
e = 0. However on reducing e gradually through its neg- 
ative values, the node with the smaller degree gets more 
preference to be connected. Therefore in the limit of 
e — > — oo the node with smaller degree is always selected 
with probability one. Even in this limit the network 
has a branched structure because if the two end nodes 
have equal degrees, any one of them gets connected with 
the new node with probability 1/2. On the other hand 
when e > 0, the node with a larger degree gets more 
weight and in the limit of e — * oo the end node with 
larger degree always gets the new connection resulting 
a star like structure. When e decreases below zero, the 
degree distribution becomes stretched exponential like: 
P(k,e < 0) ~ exp(— ak h ^) where the exponent 6(e) has 
a continuous variation with e, expected to reach 1 as 
e — > — oo and when e — > 0. Our numerical estimates 
for 6(e) for different values of e fit very well to the form 
6(e) = a (— <l) v — b Q where the constants are estimated to 
be a a w 1.19, v ps 0.14 and 6 0.85. 

To summarize, we argued that in the growing Internet 
as well as in the Airport network it is more likely that 
the new nodes get their connections in the local neigh- 
borhood. Indeed a spatial scale-free network is grown 
using the criterion of local selection rule. This network 
shows a non-linear dependence of the nodal strength on 
the degree. We conjectured that irrespective of whatever 
the dynamics may be, the non-linearity is the result of 
a power law decay of the link weight with time. When 
the same network is constructed without the underlying 
Euclidean space it gives a very efficient algorithm to gen- 
erate the Barabasi- Albert network. 
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